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ABSTRACT 



Context. Globular clusters with their large populations of millisecond pulsars (MSPs) are believed to be potential emitters of high-energy gamma- 
ray emission. The observation of this emission provides a powerful tool to assess the millisecond pulsar population of a cluster, is essential for 
understanding the importance of binary systems for the evolution of globular clusters, and provides complementary insights into magnetospheric 
emission processes. 

£T) ■ Aims. Our goal is to constrain the millisecond pulsar populations in globular clusters from analysis of gamma-ray observations. 
CD Methods. We use 546 days of continuous sky- survey observations obtained with the Large Area Telescope aboard the Fermi Gamma-ray Space 
Telescope to study the gamma-ray emission towards 13 globular clusters. 

Results. Steady point-like high-energy gamma-ray emission has been significantly detected towards 8 globular clusters. Five of them (47 Tucanae, 
r* ■ Omega Cen, NGC 6388, Terzan 5, and M 28) show hard spectral power indices (0.7 < T < 1.4) and clear evidence for an exponential cut-off 
^ t 1 in the range 1.0 — 2.6 GeV, which is the characteristic signature of magnetospheric emission from MSPs. Three of them (M 62, NGC 6440 and 
NGC 6652) also show hard spectral indices (1.0 < T < 1.7), however the presence of an exponential cut-off can not be unambiguously established. 
Three of them (Omega Cen, NGC 6388, NGC 6652) have no known radio or X-ray MSPs yet still exhibit MSP spectral properties. From the 
■ observed gamma-ray luminosities, we estimate the total number of MSPs that is expected to be present in these globular clusters. We show that 
our estimates of the MSP population correlate with the stellar encounter rate and we estimate 2600 - 4700 MSPs in Galactic globular clusters, 
commensurate with previous estimates. 

Conclusions. The observation of high-energy gamma-ray emission from globular clusters thus provides a reliable independent method to assess 
their millisecond pulsar populations. 

Key words. Pulsars: general - Globular clusters: NGC 104, NGC 5139, NGC 6266, NGC 6388, Terzan 5, NGC 6440, NGC 6441, NGC 6541, 
NGC 6624, NGC 6626, NGC 6652, NGC 6752, NGC 7078 - Gamma rays: observations 
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1. Introduction 

With their typical ages of ~ 10 10 years, globular clusters form 
the most ancient constituents of our Milky Way Galaxy. These 
gravitationally bound concentrations of ten thousand to one mil- 



lion stars are surprisingly stable against collapse which implies 
some source of internal energy that balances gravitation. The po- 
tential energy of binary systems is a plausible source of this in- 
ternal energy, tapped by close stellar encounters that harden the 
orbits of the systems dHut et al .11 19921) . Indeed, globular clusters 
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contain considerably more close binary systems per unit mass 
than the Galactic disk w hich event ually show up as rich popula- 
tions of X-ray binaries dClarkl 1 1 9751) . This scenario is strength- 
ened by the observation that the number of low-mass X-ray bi- 
nary systems containing neutro n stars is direc tly correlated with 
the stellar encounter rate (iGendre et al .112003 ah . Another conse- 
quence of this scenario is the pr esence of many millisecond pul- 
sarfl (hereafter MSPs; see e.g. ICamilo & Rasiol 120051; iRansoml 
2008), also known as 'recycled' pulsars, i.e. pulsars that were 
spun-up to millisecond period s by mass-accretion from a binary 
companion (lAlpar et al.ll 1992b . 

Observations with the Large Area Telescope (LAT) onboard 
the Fermi Gamma-ray Space Telescope ha ve confirmed MSPs 
as gamma-ray sources (Abd o et aDl2009allbh . The spectral en- 
ergy distribution of millisecond pulsars is characterised by hard 
(1.0 < r < 2.0) power law spe ctra with exponent ial cut-offs 
in the 1 - 3 GeV e nergy range (I Abdo et al.ll2009bh . Recently, 
lAbdo et all (l2009cl) presented the first detection of a globular 
cluster (GC) in the gamma-ray domain. This GC, 47 Tuc, has a 
spectral energy distribution best described by a photon index of 
1.3+0.3 with a cut-off energy of GeV (lAbdo et al Il2009cl) 

typical of the other MSPs detected to date (lAbdo et al.l l2009b). 
Further, 47 Tuc contains at least 23 MSPs, known from radio and 
X-ray observations. The lack of variability over days to months 
is consistent with MSP emission. In addition, folding the data 
on known ephemerides from the 47 Tuc pulsars reveals no sig- 
nificant detections, thus it appears that the gamma-ray emission 
is not due to a single MSP but rather attributable to an entire 
population of MSPs in this globular cluster. Using the observed, 
average efficiency of converting spin down energy into the ob- 
served gamma-ra y luminosity, cons traints can be placed on the 
MSP population dAbdo et al.H2QQ9cb . 

As the number of neutron star X-ray binaries are correlated 
with encounter rate and MSPs are the progeny of these systems, 
it would follow that the number of MSPs per globular cluster 
scales in a similar way. It is difficult to test such a correlation 
using radio and X-ray observations as the former are affected 
by dispersion and scattering by the turbulent ionized interstellar 
medium, in particular for clusters near the Galactic bulge, while 
the latter are affected by interstellar absorption rendering the de- 
tection difficult due to the low count rates observed. 

The gamma-ray domain is not affected by interstellar absorp- 
tion and there is also the added advantage that the gamma-ray 
beam s may be wider than the radio/X-ray beams (e.g. lAbdo et al.l 
l2010al) . which would permit more MSPs to be detected in the 
gamma-ray domain than those at lower energies, thus making 
gamma-ray observations ideal for testing such a correlation. 

In this paper we consider gamma-ray sources that are spa- 
tially consistent with GCs and that show the spectral character- 
istics of MSPs, i.e. that have hard power law spectra with ex- 
ponential cut-offs in the few GeV regime, that are steady, and 
that are point-like. We analyse Fermi LAT data for 13 globular 
clusters (see Table [B and include in our list the 8 globular clus- 
ters that have been formall y associated with so urces in the first 
year Fermi LAT catalogue (lAbdo et al. 2010b, hereafter named 
1FGL sources). We add two further globular clusters that lie spa- 
tially close to 1FGL sources (Omega Cen and NGC 6624), and 
include also NGC 6441 due to the high stellar collision rate tha t 
is believed to favour the formation of MSPs (iFreire et al.ll20 08), 
NGC 6752 due to its relative proximity of 4 kpc (ID' Amico et al.l 
2002|, and M 15 due to its relatively large population of known 
MSPs (lAndersonlll993h . 



1 See http://www.naic.edu/~pfreire/GCpsr.html for an updated list. 



Table 1. Globular clusters analysed in this work. 



Name 


Other name 


MbrS 


Reason for inclusion 


47 Tucanae 


NGC 104 


23 


1FGL J0023.9-7204 


Omega Cen 


NGC 5139 





close to 1FGLJ1328.2 


M62 


NGC 6266 


6 


1FGLJ1701. 1-3005 


NGC 6388 







1FGL J1735. 9-4438 


Terzan 5 




33 


1FGLJ1747.9-2448 


NGC 6440 




6 


1FGLJ1748.7-2020 


NGC 6441 




4 


high collision rate 


NGC 6541 







1FGLJ1807.6-4341 


NGC 6624 




4 


close to 1FGLJ1823.4 


M28 


NGC 6626 


12 


1FGLJ1 824.5-2449 


NGC 6652 







1FGLJ1835.3-3255 


NGC 6752 




5 


5 MSPs, nearby 


M 15 


NGC 7078 


8 


8 MSPs 



Notes. The known number of MSPs (column 3) has been taken from 
http ://w w w. naic . edu/~pf reire/GCp sr. html 



2. Observations 

2.1. Data preparation 

The data used in this work have been acquired by the LAT 
telescope aboard the Fermi Gamma-ray Space Telescope dur- 
ing continuous regular sky survey covering the period August 
8th 2008 - February 12th 2010 (546 days). Events satisfying 
the standard low-b ackground event selection ('Diffuse' events; 
lAtwood et al.ll2009h and coming from zenith angles < 105° (to 
greatly reduce the contribution by Earth albedo gamma rays) 
were used. Furthermore, we selected only events where the satel- 
lite rocking angle was inferior to 40° for the first 12 months, and 
inferior to 52° for the remaining period. We further restricted the 
analysis to photon energies above 200 MeV; below this energy, 
the effective area in the 'Diffuse class' is relatively small and 
strongly dependent on energy. All analysis was performed us- 
ing the LAT Science Tools package, which is available from the 
Fermi Science Support Center, using P6_V3 post-launch instru- 
ment response functions (IRFs). These take into account pile- 
up and accidental coincidence effects in the detector subsystems 
that were not considered in the definition of the pre-launch IRFs. 

2.2. Analysis method 

For the analysis of each globular cluster we selected events 
within squared regions-of-interest (ROIs) that have been aligned 
in Galactic coordinates. Events have been binned into 25 
logarithmically- spaced energy bins covering the range 200 MeV 
- 50 GeV and into spatial pixels of 0.1° x 0.1° in size using a 
Cartesian projection. The size and location of the ROIs have 
been chosen to avoid nearby strong sources and the bright diffuse 
emission from the Galactic plane while maintaining the largest 
possible coverage of the point-spread function for the globu- 
lar clusters of interest. The ROI definitions are summarised in 
Tabled 

The observed events in each of the ROIs have been modelled 
using components for the gamma-ray and instrumental back- 
grounds and for known sources in the field. The backgrounds 
are a combination of extragalactic and Galactic diffuse emissions 
and some residual instrumental background. The Galactic com- 
ponent has been modelled using the LAT standard diffuse back- 
ground model gll_iem_v®2 .fit for which we kept the normal- 
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Table 2. Definition of ROIs used for analysis. 



Name 


<^J2000 


<5j2000 


Size 


47 Tucanae 


" 

00 h 24 m 01.7 s 


-72°04 / 42.9 // 


10° x 10° 


Omega Cen 


13 h 26 m 45.0 s 


-47°28'37.0" 


10° x 10° 


M 62 


17 h 01 m 10.3 s 


-30°05'00.6" 


6° x 6° 


NGC 6388 


17 h 36 m 17.0 s 


-44°44'05.6" 


6° x 6° 


Terzan 5 


17 h 48 m 04.9 s 


-24°46'48.0" 


6° x 6° 


NGC 6440 


17 h 48 m 52.0 s 


-20°21'34.0" 


6° x6° 


NGC 6441 


17 h 50 m 12.0 s 


-37°03'04.0" 


6° x6° 


NGC 6541 


18 h 08 m 02.2 s 


-43°42'19.0" 


6° x6° 


NGC 6624 


18 h 23 m 40.0 s 


-30°21'38.0" 


6° x6° 


M28 


18 h 23 m 46.0 s 


-24°55'19.0" 


6° x6° 


NGC 6652 


18 h 35 m 45.0 s 


-32°59'25.0" 


6° x6° 


NGC 6752 


19 h 10 m 51.0 s 


-59°58'54.0" 


10° x 10° 


M 15 


21 h 29 m 58.0 s 


+ 12° 10' 00.4" 


10° x 10° 



Notes. Columns are (1) source name, (2) Right Ascension (J2000) of 
ROI centre, (3) Declination (J2000) of ROI centre, and (4) squared size 
ofROI. 



isation as a free parameter. The extragalactic and residual in- 
strumental backgrounds have been combined into a single com- 
ponent which has been taken as being isotropic. The spectrum 
of this component has been modelled using the tabulated model 
isotropic_iem_v®2 . txt and the normalisation has been kept 
as a free parameter0 In addition to the diffuse components, we 
included all point sources (except for the source that corresponds 
to the globular clus ter of interest) from the first-year catalogue 
(lAbdo et al.ll2010bh that were detected with TS > 25 for a given 
ROI in the background model. As usual, the Test Statistic TS 
is defined as twice the difference between the log-likelihood 
with (Xi) and without (Xo) the source, i.e. TS = 2(Xi - £o) 
(TS = 25 corresp onds to a detection significance of ~ 4cr; 
Abdo et alJ (1201 Obi) ). In general, the spectral distributions of the 
point sources have been fitted using power law models for which 
we kept the total flux and the spectral index as free parameters. 
Only if a point source has been identified as a gamma-ray pulsar, 
we fitted its spectral distribution using an exponentially cut-off 
power law for which we kept the normalisation, the spectral in- 
dex and the cut-off energy as free parameters. 

As a first analysis step, we localised the gamma-ray emis- 
sion by fitting a point source on top of the background model 
for a grid of test positions using the binned maximum likelihood 
fitting procedure gtlike. For each ROI, the grid has been cen- 
tred on the nominal position of the globular clustefl and was 
comprised of 20 x 20 positions with a spacing of 3'. We mod- 
elled the spectrum of the point source using a simple power law 
for which we kept the integrated flux and the spectral index as 
free parameters. The TS values obtained for the point source at 
the grid positions were then interpolated using a minimum cur- 
vature surface method to approximate the TS surface near the 
globular cluster. This surface was then used to localise TS max , 
the position (arj2ooo, feooo) at which TS is maximised. The 95% 
confidence region for the true source location was derived from 
the surface contour that fulfils TS max - 6.0. For simplicity, we 
approximated this contour by a circle with radius • 

As next step, we determined the spatial extent of the gamma- 
ray emission by replacing the point source in the grid search 



2 The models can be downloaded from 

http://fermi.gsfc.nasa.gov/ssc/data/access/lat/BackgroundModels.html 
3 Taken from SIMBAD. 



by 2D Gaussian-shaped intensity profiles and by repeating the 
search for various Gaussian extents cr ext . Typically, we per- 
formed the grid search for <x ext = 0.1°, 0.2°, and 0.3°, but even- 
tually we also performed runs for <x ext = 0.05° and 0.4° to bet- 
ter constrain the source extent. The dependence of TS max on 
(Text was then determined by adjusting a parabola to the pairs 
(TS max , (Text) of values, which in all cases gave a satisfactory fit 
near the maximum. Using this parabola we then determined the 
maximum TS max and the corresponding cr ext as the maximum 
likelihood estimate of the source extent. All the LAT sources 
were consistent with being point sources, and we determined the 
2<r upper extension limit from the parabola by searching the cr ext 
for which the TS decreased by 4 from its maximum value. If the 
optimum <r ex t was formally > 0, we added this (generally small) 
offset to the upper limit. 

Once we established that the gamma-ray source is spatially 
consistent with the globular cluster and that the source is not ex- 
tended, we determined the spectral characteristics of the source 
by fitting a point source at the nominal position of the globular 
cluster on top of the background model to the data. We again 
used the binned maximum likelihood fitting procedure gtlike 
for this purpose. For the spectral model of the source we used 
an exponentially cut-off power law N(E) = NqE~ y e~ E/Ec , with 
normalisation No, spectral index T, and cut-off energy E c . We es- 
timated the significance s c of the cut-off by repeating the fit us- 
ing a simple power law model, and by using s c = ^/TS C - TS p , 
where TS C and TS p are the TS values obtained using the expo- 
nentially cut-off power law and the simple power law, respec- 
tively. 

We determined 68% confidence intervals for all spectral pa- 
rameters by successively varying each parameter while optimis- 
ing the others until TS = TS max - 1.0 is fulfilled. The result- 
ing statistical parameter uncertainties are to first order consis- 
tent with the values that we obtained using gtlike, which are 
estimated from the covariance matrix of the parameters at max- 
imum likelihood. While the latter uncertainties are by definition 
symmetric about the maximum likelihood parameter values, our 
method allows for the determination of asymmetric confidence 
intervals, which are particularily relevant for describing the sta- 
tistical uncertainty in the measurements of E c . 

In addition to the statistical uncertainties we also estimated 
systematic uncertainties for all spectral parameters that arise 
from uncertainties in the precise knowledge of the effective area 
of the LAT telescope and of its energy dependence. We do this 
by repeating the gtlike fits using modified instrument response 
functions that bracket the uncertainties in the effective area. For 
each parameter we quote as systematic uncertainty the maximum 
parameter offset that is obtained with respect to the analysis done 
with the P6_V3 IRF. 

We also searched for possible flux variability of the sources 
associated with the globular cluster by fitting our source model 
to the data on a monthly basis. For this analysis we fixed the 
power law spectral indices of all sources and the normalisation 
of the isotropic diffuse component to the values we obtained for 
the full dataset. The normalisation of the Galactic diffuse compo- 
nent was kept as a free parameter due to the fact that this com- 
ponent was found to vary significantly, in particular for ROIs 
that are close to the Galactic plane. From this analysis we ob- 
tained lightcurves for all globular clusters that we fitted to a con- 
stant flux model using a x 2 minimisation procedure. Since we 
fit lightcurves for 17 months of data, we have 16 degrees-of- 
freedom in this minimisation, and the source can be considered 
as not being constant at the 68% or 95% confidence levels if 
*month > 18.1 or > 26.3, respectively. 
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Finally, we derived the spectral energy distributions (SEDs) 
of the globular clusters by dividing the energy range from 200 
MeV to 20 GeV into 8 logarithmically- spaced energy bins, and 
by fitting the fluxes or normalisations of all model components 
independently for each bin. For all sources, we assumed sim- 
ple power law models with a fixed spectral index of Y = 2.1 
within each energy bin. Since the energy bins chosen are rel- 
atively narrow, the precise value of the slope has little impact 
on the resulting spectra. The normalisations of the diffuse com- 
ponents (Galactic diffuse and isotropic backgrounds) are fitted 
independently for each energy bin, which makes this analysis 
less sensitive to systematic uncertainties in the spectral energy 
distribution of these model components. 

We also note that we searched for pulsations from MSPs 
in globular clusters and, to date, no significant pulsations have 
been reported from any of the pulsars that has been searched. 
However, we only have accurate ephemerides (derived from 
times of arrival measured with the Parkes, Nancay, and Jodrell 
Bank radio telescopes) for 15 globular cluster MSPs, limiting 
our search to only ~ 10% of the known globular cluster pul- 
sar population. This illustrates the pressing need for more con- 
temporaneous timing solutions in order to reveal the gamma-ray 
pulsations from globular cluster MSPs. 

2.3. Results 

For 11 of the 13 ROIs analysed, a gamma-ray source has been 
found in the Fermi LAT data that is spatially consistent with a 
globular cluster. For two ROIs (NGC 6441 and NGC 6624) we 
were unable to find a convincing gamma-ray counterpart; how- 
ever, gamma-ray sources were detected close by to both globular 
clusters. For NGC 6441 we found evidence for a faint gamma- 
ray source (TS = 29) that we localised at #j2ooo = 17 h 50.5 m and 
feooo = -36°52 r with r 95 = 7.5', yet with an offset of 11.7', 
the globular cluster is located just outside the 99% error contour . 
NGC 6624 is close to 1FGL J1823.4-3009 (lAbdo et al.ll20l0bh 
which we detected with TS = 43, and which we localised at 
tfj20oo = 18 h 23.5 m and £ J2 ooo = -30°13' with r 95 = 6.3'. With 
an offset of 8.6', NGC 6624 is also located just outside the 99% 
error contour of the gamma-ray source. Formally, we can not 
establish an association of NGC 6441 and NGC 6624 with the 
gamma-ray sources, and hence we excluded both objects from 
further analysis. Both gamma-ray sources, however, are rather 
faint, and more data are needed before an association of these 
sources with globular clusters definitely can be excluded. 

Three of the 11 gamma-ray sources that are spatially con- 
sistent with a globular cluster have TS < 25 (NGC 6541, 
NGC 6752, and M 15), hence we cannot claim a significant 
detection for these cases. NGC 6541 has been a ssociated in 
the firs t year catalogue with 1FGL J1807.6-4341 (lAbdo et al.1 
1201 Obi) , yet with TS = 25.3 this source was just barely above the 
threshold of TS = 25 used for compilation of the catalogue. In 
our dedicated analysis we find a lower detection significance of 
TS = 12.0 for the source, which can be explained by differences 
in the analysis method and data selection with respect to the cat- 
alogue analysis. In particular, the catalogue analysis determined 
TS using an unbinned maximum likelihood analysis which has 
recently been recognized to provide larger TS values in several 
low-latitude fields than the binned analysis 

This leaves us with 8 significant globular cluster candidates 
for which we show smoothed raw and residual counts maps in 
Fig.Q] For each of these clusters we determined fluxes and spec- 
tral parameters using exponentially cut-off power laws. For the 
3 non- significant globular clusters we derived formal 2cr upper 



flux limits by fitting a point source at the nominal cluster posi- 
tion to the data for which we increased the flux, starting from 
its maximum likelihood value, until the TS decreased by 4 with 
respect to the maximum. As a spectral model we used an expo- 
nentially cut-off power law for which we fixed the spectral index 
and the cut-off energy to the weighted average values of T = 1.4 
and E c - 1.6 GeV that were obtained for the 8 significantly de- 
tected globular cluster candidates. We added to these statistical 
upper limits the systematic uncertainty that has been obtained 
using the bracketing IRFs (cf. Section [2T2b and that accounts for 
uncertainties in the effective area of the LAT. The fluxes and 
spectral characteristics as well as the maximum likelihood lo- 
cations of the gamma-ray sources are summarised in Table [3] 
Quoted uncertainties are statistical and systematic. 

Table [3] quotes also the variability indicator A^ onth , which is 
inferior to 26.3 for all globular cluster candidates, indicating that 
none of the gamma-ray sources exhibit significant time variabil- 
ity. All 8 sources thus match the expected property of globular 
clusters for being non- variable sources of gamma rays. Table [3] 
also shows that the 8 significant sources exhibit the expected 
spectral characteristics of globular clusters: a flat (r < 2) power 
law spectral index with an exponential cut-off in the few GeV 
domain. For two of the sources (M 62 and NGC 6440), how- 
ever, the statistical significance of the cut-off is s c < 3<x. More 
data are clearly needed here before definite conclusions about 
the reality of the exponential cut-off can be drawn. Finally, we 
show in Fig. [2] the spectral energy distributions of the 8 signif- 
icant globular cluster candidates. In the following we comment 
on the analysis results for each of these globular clusters. 

2.4. Comments on individual clusters 
2.4 A . 47 Tucanae 

47 Tuc was the first globular cluster to b e associated with a 
gamma-ray source detected by Fermi LAT (lAbdo et al.i r2009c). 
The maximum likelihood position of the gamma-ray source is 
only 1.7' offset from the cluster core, which is well inside the 
95% confidence error radius of 3.3'. The spatial and spectral 
characteristics that we derive from the present datas et are com- 
patible with those published in our previous work (lAbdo et al.l 
2009c), yet the ~ 2.8 times longer exposure now allows us 
to constrain the spectral source parameters considerably bet- 
ter. In particular, the uncertainty on the energy flux of 47 Tuc 
has been reduced by a factor ~ 2 to less than 10%, improving 
the constraint on the total gamma-ray luminosity of the clus- 
ter, and hence on the expected population of millisecond pulsars 
(cf. Section [3]). 

2.4.2. Omega Cen 

Omega Cen is one of the globular clusters in our list for 
which no MSP has so far been detected. Omega Cen is situ- 
ated near the southern giant radio-lobe of the nearby Cen A 
radio galaxy which is a lso an extended source of GeV gamma 
rays (lAbdo et al.l 1201 Pel) . To account for emission from the ra- 
dio lobes in our analysis we included a s patial template based 
on the 22 GHz WMAP microwave map (IHinshaw etall 12009b 
in our background model. Formally, the position of the gamma- 
ray source near Omega Cen (1FGL J1328. 2-4729) is inconsis- 
tent with the location of the globular cluster. It turned out, how- 
ever, that 1FGL Jl 328.2-4729 is a superposition of two distinct 
sources that can be separated thanks to their different spectral 
characteristics. While the first source shows a cut-off in the few 
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Fig. 1. Gaussian kernel smoothed (<x = 0.1°) counts maps of the ROIs for the 8 significant globular cluster candidates. The insets 
show adaptively smoothed background subtracted counts maps for a 3° x 3° large region around the globular cluster. Circles indicate 
the nominal locations of the globular clusters. 



GeV domain, the second source is rather hard (T = 1.6 + 0.2), 
with no cut-off and can be clearly isolated from the first one by 
selecting only events > 10 GeV. This results in a best fitting loca- 
tion of a J2 ooo = 13 h 28.7 m and £ J2 ooo = -47°30' with a 95% error 
radius of = 4.0'. Alternatively, we fixed a point source at 
the location of Omega Cen and fitted the location of the residual 
counts to anooo = 13 h 28.8 m and S moo = -47°29' with a 95% 
error radius of = 3.1'. Both locations are spatially consistent 
with the radio source SUMSS Jl 32840-472748, which has for- 
mally been associated to 1FGL Jl 328.2-472 9 in the Fermi LAT 
first year catalogue of active galactic nuclei (lAbdo et al.ll2010dl) . 
The gamma-ray source is also consistent with the radio source 



PMN J1328-4728 and the X-ray sources 1WGA J1328.6-4727 
and 1RXS Jl 32846.7-472759, which all are possibly just dif- 
ferent designators for the same source detected by other instru- 
ments. 

We thus include this source in our background model and 
re-determine the localisation of the source that is only visible 
below a few GeV. This results in a maximum likelihood position 
of ttnooo = 13 h 26.5 m and feooo = -47°29' with a 95% error 
radius of rgs - 7.5'. This position is 3.2' offset from the centre 
of Omega Cen, and thus the gamma-ray source is spatially fully 
consistent with the location of the globular cluster. 
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Table 3. Gamma-ray characteristics of globular clusters. 



Name 


<*J2000 


<5j2000 


r 9 5 


°~ ext 


TS 


% month 


Photon flux 


Energy flux 


r 


E c 


S c 


47 Tucanae 


00 h 


23.8 m 


-72°04' 


3.3' 


< 


4.8' 


603.3 


9.6 


o Q+0.6+0.4 
-0.5-0.3 


9 c +0.2+0.2 
-0.2-0.2 


i 4+0.2+0.2 
-0.2-0.2 


9 9+O.8+O.3 
-0.5-0.2 


5.6 


Omega Cen 


13 h 


26.5 m 


-47°29' 


7.5' 


< 


8.4' 


50.0 


14.6 


a q+0.5+0.3 
U ' ^-0.4-0.2 


i A+0.2+0.1 
1,W -0.2-0.1 


c\ 7+0.7+0.4 
u ' 7 -0.6-0.4 


i 9+0.7+0.2 
1 '^-0.4-0.2 


4.0 


M62 


17 h 


01. l m 


-30°08' 


4.4' 


< 


7.2' 


107.9 


16.0 


9 7+1.0+1.9 
z " ' -0.9-0.8 


9 i +0.3+0.5 
z " 1 -0.3-0.1 


i 7+0.3+0.4 
i ' 7 -0.3-0.5 


A 4+3.8+17.7 
^ -1.8-1.8 


2.5 


NGC 6388 


17 h 


35.9 m 


-44°41' 


5.7' 


< 


9.0' 


86.6 


13.8 


i /r+1.0+2.0 
i ' u -0.6-0.6 


i ^+0.3+0.6 
1,u -0.3-0.2 


i 1+0.7+0.8 
x ' -0.5-0.8 


i + 1.2+1.8 
i '°-0.7-0.6 


3.3 


Terzan 5 


17 h 


47.9 m 






< 


9.0' 




ZD. J 


7 zr + 1.7+3.4 
' -1.5-2.2 


7 i +0.6+1.0 
7 ' -0.5-0.5 


i 4+0.2+0.4 
1 -0.2-0.3 


9 6+0-7+1.2 
^•"-0.5-0.7 


7 1 
/ . 1 


NGC 6440 


17 h 


48.8 m 


-20°21' 


5.2' 


< 


8.4' 


65.7 


5.9 


9 q+2.7+4.4 
^•^-1.3-1.1 


9 9+0.9+1.2 
^-0.5-0.2 


i 6+0-5+0.6 
i ' u -0.5-0.8 


3 1+ 3.3+oo 
J -1.4-1.1 


1.4 


M28 


18 h 


24.4 m 


-24°51' 


8.0' 


< 15.6' 


77.9 


20.6 


9 £+1.3+2.2 
^ -1.0-0.9 


9 a+0.4+0.6 
^-0.3-0.3 


i i +0.7+0.6 
A ' -0.5-0.7 


i rj+0.6+0.4 
1 -0.3-0.2 


4.3 


NGC 6652 


18 h 


35.7 m 


-33°01' 


7.5' 


< 


9.6' 


54.8 


9.8 


r> 7+0.5+0.2 
u ' 7 -0.3-0.1 


n 0+0.2+0.1 


i n+0.6+0.3 
1,w -0.5-0.3 


i O + 1.2+0.4 
i '°-0.6-0.3 


3.2 


NGC 6541 
NGC 6752 
M 15 


18 h 07.9 m 
19 h 10.3 m 
21 h 29.4 m 


-43°41' 
-59°56' 
+ 12W 


20.1' 
6.3' 
6.9'* 






12.0 
13.7 
5.4 




< 1.1 

< 0.7 

< 0.6 


<0.8 
<0.5 
<0.5 


(1.4) 
(1.4) 
(1.4) 


(1.6) 
(1.6) 
(1.6) 





Notes. Columns are (1) source name, (2) fitted Right Ascension (J2000) of LAT source, (3) fitted Declination (J2000) of LAT source, (4) 95% 
confidence level error radius (*for M 15 the 68% confidence level error radius is quoted because the small significance of the detection does not 
allow to define a meaningful 95% confidence level error radius), (5) 2cr upper limit on source extent, defined as cr of a 2D Gaussian intensity 
profile, (6) value of Test Statistic, (7)^ onth as a measure of source variability, (8) integrated > 100 MeV photon flux in units of 10~ 8 ph cm -2 s" 1 , 
(9) integrated > 100 MeV energy flux in units of 10 -11 erg cm -2 s _1 , (10) spectral index T of exponential cut-off power law, assuming N(E) oc 
E~ r e~ E/Ec , (11) cut-off energy E c in GeV, and (12) significance of exponential cut-off in units of Gaussian cr. 



2.4.3. M 62 



M 62 has formally been associated with 1FGL J1701. 1-3005 
(lAbdo et al.ll2010bl) and also our maximum likelihood position 
is only 1.6' offset from the core of the globular cluster, which is 
thus well within = 4.4'. The gamma-ray source shows only 
marginal (s c = 2.5) evidence for an exponential cut-off and ex- 
hibits the largest formal cut-off energy (E c = 4.4+-* '^} 7 8 7 GeV) 
of all sources. The SED (cf. Fig. [2]) shows no clear indication 
for a distinct cut-off, and if there were not a lack of flux be- 
low ~ 600 MeV, the SED would probably also be compatible 
with a soft power law. M 62 is located in a region of intense 
Galactic diffuse emission (angular distance from the Galactic 
centre < 10°; cf. Fig.Q]) which makes any proper extraction of 
spectral parameters difficult. In particular, we recognised that 
the source spectrum is rather sensitive to the exact choice of 
the ROI, which is readily explained by systematic uncertainties 
in the model of Galactic diffuse emission that we used for the 
analysis. We thus qualify the gamma-ray source associated with 
M 62 only as possible globular cluster candidate (in contrast to 
being a plausible candidate), and await the accumulation of more 
data before drawing any definite conclusions. 



2.4.4. NGC 6388 



NGC 6388 has formally been associated with 
1FGL J1735.9-4438 (lAbdo et al] l2010bh . Our position is 
consistent with that of 1FGL J1735.9-4438 and NGC 6388 is 
found right on the 95% confidence contour of the gamma-ray 
source, at an angular separation of 5.7'. The SED shows a clear 
cut-off at E c = GeV whic h is typical of the c ut-off 

energies that are observed for MSPs (lAbdo et al . 2009b). The 
Fermi LAT source thus qualifies as a plausible candidate for a 
globular cluster. We note that no MSP has so far been detected 
in this GC. 



2.4.5. Terzan 5 

Terzan 5 has formally been associated with 1FGL J1747. 9-2448 
(lAbdo et al.ll2010bl) and our maximum likelihood position is off- 
set by 2.4' from the cluster centre, which is inferior to r^ = 2.9' 
fo r this source. 1FGL J1747. 9-2448 has recently been studied 
by Kong et al. (2010). Using an exponentially cut-off power law 
model, they obtained a spectral index of T = 1.9 ± 0.2, a cut- 
off energy of E c = 3.8 ± 1.2 GeV, and 0.1-10 GeV pho- 
ton and energy fluxes of 2 x 10" 7 ph cm -2 s" 1 and 1.2 x 10~ 10 
erg cm -2 s" 1 , respectively. Our analysis confirms the presence of 
the exponential cut-off at a significance of 7. 1<t, yet with slightly 
different spectral parameters. In particular, our fluxes are lower 
by about a factor of 2, and our SED (our Fig. [2]) more closely 
follows an exponentially c ut-off model rather than the SED de- 
rived by lKong et al] (120101) (their Fig. 2). Possibly, the choice of 
a rather large ROI with a radius of 15° for their analysis intro- 
duced some systematic uncertainties due to the inaccurate sub- 
traction of Galactic diffuse emission in this complex region near 
the Galactic centre. We used a rectangular ROI of 6° x 6° for our 
analysis which considerably reduces the impact of the Galactic 
diffuse model on the analysis, and which allows for a more reli- 
able determination of the spectral parameters of the source. The 
spectral parameters we obtain for Terzan 5 are indeed very close 
to those observed for 47 Tuc, making the gamma-ray source a 
very plausible candidate for being a globular cluster. 

2.4.6. NGC 6440 

NGC 6440 has formally been associated with 
1FGL J1748.7-2020 (lAbdo et al j 1201 Obh . and our fitted 
position is offset by 1.3' from the cluster core, well within 
= 5.2 r . The SED, however, shows no convincing evidence 
for an exponential cut-off, and also the spectral analysis results 
in only s c = 1 .4<r, too small to claim the detection of a cut-off. 
The spectral analysis may have suffered from the proximity 
of the Galactic plane and the related uncertainties due to the 
subtraction of the diffuse emission, as well as th e proximity 
of a bright gamma-ray pulsar (PSR J1741-2054: lAbdo et al] 
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Fig. 2. Spectral energy distributions of the 8 significant globular cluster candidates. Lines indicate the exponentially cut-off power 
law models that have been fitted to the data using the binned maximum likelihood fitting procedure, and for which the spectral 
parameters are quoted in Table [3] The dotted part of the lines indicate the extrapolation of the exponentially cut-off power law 
model towards energies below 200 MeV. 



(2009d)). In any case, the gamma-ray source associated with emission. Pending the acquisition of these data we thus qualify 
NGC 6440 is a relatively faint source, and more data are the source only as a possible candidate for being a globular 
required to allow for an accurate spectral characterisation of the cluster. 
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2.4.7. M 28 

M 28 has been for mally associated with 1FGL Jl 824.5-2449 
(lAbdo etal.1 1201 Obi) , and our source position is offset by only 
1.6' from the centre of the cluster, well within the 95% error 
radius of 8.0'. Figure [T] suggests a somewhat elongated shape for 
the gamma-ray source that may indicate a superposition of two 
objects, yet the source is formally consistent with a point source, 
although with the largest of all upper limits of <x ex t < 15.6'. The 
SED of M 28 shows a clear exponential cut-off, which is detected 
at a significance level of 4.3<x, the third largest significance after 
Terzan 5 and 47 Tuc. The gamma-ray source associated to M 28 
thus qualifies as a plausible globular c luster candidate. 

Note that iPellizzoni et al.1 (120091) have reported the detec- 
tion (at 4.2cr significance) of pulsations from PSR J1824-2452A 
with the AGILE satellite during a 5 -day interval in August 
2007, while no pulsations have been seen for other periods 
of comparable duration. PSR J1824-2452A is one of the 12 
known MSPs in M 28 and ranks among the y oungest and 
most energetic millisecond pulsars that are known (Knight et al.l 
20() 6j). The >!() ( ) MeV flux of (18 ± 5) x 10" 8 ph cm" 2 s" 1 that 
Pelli zzoni et al.l d2009) quote for the 5 -day period is significantly 
larger than the flux of (2.6 ± 1.3) X 10" 8 ph cm -2 s" 1 that we 
measured from the cluster using our dataset. Searches for pul- 
sations from PSR J1824-2452A in the Fermi LAT data have so 
far not revealed any significant detection. Furthermore, we did 
not find any evidence for significant flux variations from M 28 
on time- scales from 1 week to 1 month. Typical flux measure- 
ment uncertainties for 1 week, however, amount to ±8 x 10" 8 
ph cm -2 s" 1 (lcr), hence flares from M 28 of the duration and 
amplitude quoted by AGILE formally cannot be excluded by our 
data. 

2.4.8. NGC 6652 

NGC 6652 has been formally associated with 
1FGL J1835.3-3255 (lAbdo et al.l l2010bl) . and our source 
position is offset by 1.7' from the cluster, well below - 7.5'. 
The SED does not reveal a convincing exponential cut-off 
although the spectral analysis suggests s c = 3.2<x, mainly due 
to the lack of flux below ~ 600 MeV. The gamma-ray source 
associated with NGC 6652 is a rather faint source, and also 
here, more data are required to better constrain the SED before 
definite conclusions should be drawn. Pending the acquisition 
of these data we thus qualify the source only as a possible 
candidate for being a globular cluster. Establishing a firm 
exponential cut-off is particularily interesting in this case since 
so far no MSP has been detected in NGC 6652. 



3. Discussion 

3.1. Omega Cen, NGC 6388 and NGC 6652 

Three of the globular clusters for which we detected gamma-ray 
counterparts, Omega Cen, NGC 6388 and NGC 6652, deserve 
attention, since so far no MSPs have been detected in these ob- 
jects. Present pulsar searches in the radio and X-ray domains 
are severely limited by the instrument sensitivity, and only in a 
few examples (e.g. 47 Tuc and M 15) has th e end of the pul- 
sar lu minosity function possibly been reached (ICamilo & Rasiol 
2005]). In some cases, this may simply be because pulsars are 
generally intrinsically weak objects, and globular clusters are 
often distant. Interstellar scattering, which broadens pulsations 
because of multipath propagation, can also be a major obsta- 



cle, especially for MSPs in clusters at low Galactic latitudes. 
Furthermore, interstellar absorption hinders detection in the X- 
ray range, in particular for clusters that are situated near the 
Galactic plane. 

NGC 6388 and NGC 6652 are indeed the most distant ob- 
jects in our list of globular clusters for which gamma-ray emis- 
sion is significantly detected (see Table©. Among the 26 globu- 
lar clusters in which MSPs have so far been detectecQ, 6 have dis- 
tances that are comparable to or larger than those of NGC 6388 
and NGC 6652. Four of these clusters (NGC 1851, M 53, M 3, 
and M 15) show considerably lower dispersion measures than 
NGC 6388 and NGC 6652, hence radio pulsars are more eas- 
ily detected in these systems despite their large distances. The 
two other clusters (NGC 6441 and NGC 6517) have dispersion 
measures that are comparable to the estimates for NGC 6388 and 
NGC 6652, yet uncerta inties in the dispersion measure estimates 
dCordes & Laziol 120021) or selectio n effects in the radi o surveys 
that are very difficult to quantify (Hessels et al. 2007) may ex- 
plain why no MSPs have so far been detected in these clusters. 

Omega Cen is substantially closer than NGC 6388 and 
NGC 6652, at a d istance that is comparable to that of 47 Tuc 
dBono et al .1120081) . However, while 47 Tuc is situated at high 
galactic latitudes with a correspondingly small dispersion mea- 
sure of ~ 25 pc cm -3 , Omega Cen is closer to the galac- 
tic pla ne and has an estim ated dispersion measure of 126 pc 
cm -3 dCordes & Laziol 12002), so MSPs searches will be more 
difficult in Omega Cen compared t o 47 Tuc. Om e ga C en has 
been searched for radio pulsars by lEdwards et al] (120011) with 
the Parkes 64 m radio telescope at 660 MHz without suc- 
cess. Current X-ray observations of this cluster with Chandra 
and XMM-Newton ( Haggard et all 120091: iGendre et al.1 l2003bl) 
are not deep enough to uniquely identify MSPs. However, 
Has sard et al. (2009) note that their Chandra observations do 
not rule out the presence of MSPs, as numerous sources occupy 
the region of the colour-magnitude diagram (CMD) where MSPs 
would be expected to l i e. The y compared their X-ray CMD with 
that of iGrindlav etal.1 (120011) for 47 Tuc and adjusted for the 
somewhat larger hydrogen column density to Omega Cen ver- 
sus 47 Tuc (and slightly different distances and exposure times). 
Eighteen objects are found in the range of brightness and colour 
expected for MSPs in Omega Cen. This range of colours is also 
typical of chromospherically active single stars and binaries, but 
many or all of these objects may be MSPs. Our estimate is also 
commensurate with the 28 (with an r ms dispersion of 6) MSPs 
pred icted via two body tidal capture (iDi Stefano & R appaport 
11992b . 



3.2. Expected MSP populations 

Under the assumption that MSPs share the same characteristics 
in the various globular clusters that we studied, we may use the 
observed gamma-ray luminosities to obtain estimates on the to- 
tal number of MSPs expected to reside in each cluster. Here 
we closely follow the formalism that has been developed by 
lAbdo et all (l2009cl) to estimate the size of the MSP population 
in 47 Tuc. The total number Nmsp of MSPs present in a globu- 
lar cluster is estimated from the isotropic gamma-ray luminosity 
of the cluster L r , the average spin-down power (E) of MSPs, 



4 See http://www.naic.edu/~pfreire/GCpsr.html 
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Table 4. Isotropic gamma-ray luminosities and expected num- 
bers of MSPs. 



Name 


d (kpc) 


L r (10 J erg s ) 


A^MSP 


47 Tucanae 


4.0 + 0.4 (1) 


4.8!i;j 




Omega Cen 


4.8 + 0.3 (2) 




19+1 


M62 


6.6 + 0.5 (3) 


10 - 9 -23 


/u -34 


NGC 6388 


11.6 + 2.0 (4) 


25.8 + !^ 

-10.6 


180 + ! 2 ° 

-100 


Terzan 5 


5.5 + 0.9 (5) 


25.7^ 


180^° 


NGC 6440 


8.5 + 0.4 (6) 


19 +131 


130 +10 ° 


M28 


5.1 +0.5 (7) 




4 o+24 
* D -21 


NGC 6652 


9.0 + 0.9 (8) 


7.8!g 




NGC 6541 


6.9±0.7 (9) 


<4.7 


<47 


NGC 6752 


4.4 + 0.1 (10) 


< 1.1 


< 11 


M 15 


10.3 + 0.4 (11) 


< 5.8 


<56 



Notes. References to the distance estimates are: mlMcL aushlin et al. 



(2006 ); (2) 



van de Yen et all (120061): (3 ) Beccari et al l (120061): (4) 
iMoretti et all (120091): (5) lOrtolani et"aLl (120071): (6) lOrtolani et all 
(I1994D: (7) IRees & Cudworthl (119911): (8) IChabover et all (12000): 



(9) "lLee & Carnevl (120 06); flO) ISaraiedini et all d2007h : (11) 
Ivan d en Bosch et al. (200g). 



For NGC 6440. lFreire et al.l (120081) estimate that it may house as 
many pulsars as Terzan 5, a pr e dictio n that is compatible with 
our observations. Hein ke etaD (2005) estimate that a total of 
-25 MSPs exist in 47 Tuc (<60 at 95% confidence), regardless 
of their radio beaming fraction, based on their deep Chandra 
X-ray observations of this cluster This is comparable to t he up- 
per limit of 30 MSPs derived bv lMcConnell et al.l (120041) from 
the radio lum i nosity distr ibution of pulsars. As n o ted ab ove, 
Haggard et al] (120091) and |DLStefano & Rappaport (1992) ob- 
tain s imilar estimate s for th e number of MSPs in Omega Cen as 
well. KulkarnTetaD (Il990h used radio observations of Galactic 
globular clusters visible from the northern hemisphere and took 
into account selection effects in the various radio surveys, and 
then after estimating the relative efficiency of pulsar production 
in the individual clusters, produced a census of the cluster pop- 
ulation of pulsars. Their estimates are surprisingly close to ours 
where they deduce 150 MSPs in Terzan 5, 86 in NGC 6440, 18 
in M 28, 21 in M 15 but only 5 in NGC 6652, the latter being the 
only cluster for which they estimate a quantity that is remark- 
ably dissimilar to our work (see Table |4]). We recall, however, 
that we classified NGC 6652 only as a possible globular cluster 
candidate because of the rather weak evidence for an exponential 
cut-off in the spectrum. The acquisition of more data is definitely 
required before any firm conclusions about the number of MSPs 
in NGC 6652 should be drawn. 



and the estimated average spin-down to gamma-ray luminosity 
conversion efficiency (n y ) using 



Ni 



MSP 



(1) 



We compute the isotropic gamma-ray luminosity using L y = 
AnS d 2 , where S is the observed energy flux and d the distance 
to the globular cluster. (E) for the pulsar populations in individ- 
ual clusters are only poorly known since measurements of E are 
significantly affected by ac celeration in the grav itational poten- 
tial of the globular cluster (Bogdanov et al. 2006). We therefore 
adopt a value of (E) = (1.8±0.7)xl0 34 erg s" 1 for all clusters as a 
best estimate that has been derived by comparing the log P distri- 
bution of Galactic field MSPs to the acceler ation corrected log P 
distribution obtained for MSPs in 47 Tuc (lAbdo et al.ll2009cl) . 
This value is close to the value of (E) = 3.0 x 10 34 erg s" 1 
that is obtained as average for Galactic field MSPs from the 
ATNF catalogue. Since intrinsically luminuous pulsars are more 
easily detected, this average is likely biased towards high val- 
ues. Limiting the average to the local (< 1 kpc) Galactic field 
MSPs in the ATNF catalogue indeed reduces the average value 
to (E) = 1.1 x 10 34 erg s" 1 , consistent with our adopted value. 
Furthermore, we assume (rj y ) = 08 as derived from observa- 
tions of nearby MSPs (lAbdo et al.ll2009cl) . The results are sum- 
marised in Table |4] The statistical and systematic uncertainties 
in the energy flux have been added in quadrature and we quote 
the resulting uncertainty from that combination for L y and Amsp- 
It is encouraging to note that the estimates of MSP popula- 
tions in the above globular clusters are very similar to theoreti- 
cal and/or observational constraints made in other wavelengths. 
In Terzan 5, 33 MSPs are known so far and many more are ex- 
pect ed based on the extra polation of the radio luminosity func- 
tion dRansom et al.l l2005). Indeed, by comparing the total radio 
luminos ity of Terzan 5 to the ra dio luminosity function of field 
pulsars. iFruchter &~Go ss ( 2000) esti mate the number o f pulsars 
in the cluster to be 60 - 200, and also Ko ng et all ([2010) estimate 
that the number of pulsars in Terzan 5 could be as large as ~ 200. 



3.3. Stellar encounters and the formation of MSPs 

It has been shown for globular clusters that the number of neu- 
tron star X-ray binaries i s correlated with the stellar encounter 
rate T e (see Section 1 and lGendre et al .112003 al) . Since MSPs are 
believed to be the progeny of these systems, we may expect a 
similar correlation between our estimated number of MSPs (or 
equivalently the observed gamma-ray luminosity) and T e . 

To test the link between the formation of MSPs and the stel- 
lar encounter rate, we compute the latter using T e = Pq' 5 ^, where 
po is the central cluster density, given in units of L pc -3 , and r c 
is the cluster core radius, given in units of pc. We used t he cen- 
tral c luster densities and angular core radii 6 C quoted in Harris 
(119961) (February 2003 revision) along with the di stance esti- 
mates listed in Table |4] using po = p^ arns d narns /d dDjorgovskil 
Il993l) . and using r c = dtan# c . We then normalised the stellar 
encounter rates so that T e = 100 for M 62 (for M 62 we obtain 
p l 5 r 2 = 6.5 x 10 6 Lq 1 ' 5 pc~ 2 - 5 ). 

Figure [3] shows evidence for a linear correlation between 
the predicted numbers of MSPs and the stellar encounter rates. 
Taking into account the uncertainties, the correlation is well fit- 
ted by Nmsp = (0.5 ± 0.2) x T e + (18 ± 9) (correlation coefficient 
0.7). A similar correlation (not shown) is obtained if we plot the 
isotropic gamma-ray luminosity versus the encounter rate, since 
A^msp and L y are linearly related by Eq. (Q]). The correlation is 
well fitted by L y = (5.5±1.9)xl0 32 xr e +(2.6±0.7)xl0 34 erg s" 1 . 

Using our correlation we predict 2600 - 4700 MSPs that 
are observable in ga mma rays in Galactic globular clusters. 
Kul karni et all (Il990h estimate a total of ~ 10,000 pulsars 
in Galactic globular clusters, subject to uncertainties in pul- 
sar beaming and binarity. It has been suggested that the rare 
electron-capture accretion induced collapse channel will have 
a low kick velocity and may thus d ominate the final cluster 
MSP population (llvanova et al. 112008b . However, since binaries 
dominate the production channels, the extrapolation from the 
observed sample to the total population based on encounter 
rate, used for the estimate of the total pulsar population by 
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Fig. 3. Predicted number of MSPs versus stellar encounter rate 
T e . Horizontal error bars indicate the uncertainties in T e due to 
the distance uncertainties given in Tableland due to uncertain- 
ties in 6 C that we estimated from the spread of values quoted in 
the recent literature. The data have been fitted by a linear relation 
Af MS p = 0.5xr e + 18. 



Kulka rni et all d 1990b . should still be relevant. On the low end, 
Hei nke et al.l (12005 ) estimated the Galactic globular cluster pul- 



sar population as 700, essentially from X-ray observations and 
the stellar encounte r rate, the latter being commensurate with 
the estimate of lWijers & va n Paradii s (1 1 99 ll) which was deduced 
from radio observations. iFruchter & Gossl (Il990l) deduced the 
total number of MSPs in the Galactic globular system to lie be- 
tween 500 and 2000, but concluded that the expected number of 
pulsars in a globular cluster depends only weakly on the stellar 
collision rate. 

It appears that our independent method of determining the 
number of MSPs in Galactic globular clusters through gamma- 
ray observations is entirely compatible with earlier estimates. 



4. Conclusions 

An analysis of Fermi LAT data from 13 globular clusters has 
revealed 8 significant, point-like and steady gamma-ray sources 
that are spatially consistent with the locations of the clusters. 
Five of them (47 Tuc, Omega Cen, NGC 6388, Terzan 5, and 
M 28) show hard spectral power law indices (0.7 < T < 1.4) and 
clear evidence for an exponential cut-off in the range 1.0 - 2.6 
GeV, which is the characteristic signature of magneto spheric 
emission from MSPs. We thus classify these 5 sources as plausi- 
ble globular cluster candidates. Three of them (M 62, NGC 6440 
and NGC 6652) also show hard spectral indices (1.0 < T < 1.7), 
however the presence of an exponential cut-off cannot unam- 
biguously be established. More data are required before definite 
conclusions can be drawn; hence we qualify these 3 sources as 
possible globular cluster candidates. 

From the 8 globular clusters that are associated with sig- 
nificant gamma-ray sources, 5 are known to harbour MSPs. In 
Omega Cen, NGC 6388 and NGC 6652, however, no MSPs have 
so far been detected, neither by radio nor by X-ray observations. 
The observation of gamma-ray signatures that are characteris- 
tic of MSPs provides strong support that these GCs indeed also 
harbour important populations of MSPs. In particular, we pre- 
dict from the observed gamma-ray luminosities that the total 
MSP populations amount to 10 - 30 (Omega Cen), 80 - 300 
(NGC 6388), and 30 - 80 (NGC 6652) in these clusters. Deep 



radio and X-ray follow-up observations may help to unveil first 
members of these populations. 

Our predicted number of MSPs shows evidence for a posi- 
tive correlation with the stellar encounter rate in a similar way 
to their progenitors, the neutron star low mass X-ray binaries. 
This correlation allows us to deduce the total number of MSPs 
in Galactic globular clusters (2600 - 4700) which lies midway 
between all previous estimates, supporting such a correlation. 
Such an estimate can be used to derive constraints on the original 
neutron star X-ray binary population, essential for understanding 
the importance of binary systems in slowing the inevitable core 
collapse of globular clusters. 
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